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Open 



Disturbances act as powerful structuring forces on ecosystems. To ask whether environmental 
microbial communities have capacity to recover after a large disturbance event, we conducted a 
whole-ecosystem manipulation, during which we imposed an intense disturbance on freshwater 
microbial communities by artificially mixing a temperate lake during peak summer thermal 
stratification. We employed environmental sensors and water chemistry analyses to evaluate the 
physical and chemical responses of the lake, and bar-coded 16S ribosomal RNA gene pyrosequen- 
cing and automated ribosomal intergenic spacer analysis (ARISA) to assess the bacterial 
community responses. The artificial mixing increased mean lake temperature from 14 to 20 C for 
seven weeks after mixing ended, and exposed the microorganisms to very different environmental 
conditions, including increased hypolimnion oxygen and increased epilimnion carbon dioxide 
concentrations. Though overall ecosystem conditions remained altered (with hypolimnion 
temperatures elevated from 6 to 20 C), bacterial communities returned to their pre-manipulation 
state as some environmental conditions, such as oxygen concentration, recovered. Recovery to 
pre-disturbance community composition and diversity was observed within 7 (epilimnion) and 11 
(hypolimnion) days after mixing. Our results suggest that some microbial communities have 
capacity to recover after a major disturbance. 
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Introduction 

Disturbances, defined as events that trigger a 
discrete change in the physical or chemical envir- 
onment that may affect the local community (Glasby 
and Underwood, 1996), have a key role in the 
structuring of ecosystems (Sousa, 1984; White and 
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Pickett, 1985). Characteristics and ecological 
consequences of disturbances vary immensely 
(Fraterrigo and Rusak, 2008), hindering advances 
in the theoretical and predictive models (Sousa, 
1984; White and Jentsch, 2001). As there is an 
inverse relationship between disturbance frequency 
and magnitude (large events seldom occur, small 
events occur frequently (Sousa, 1984; Romme et ah, 
1998; Turner et ah, 1998; Turner and Dale, 1998; 
White and Jentsch, 2001), extreme disturbances are 
inherently difficult to study in nature. Nonetheless, 
observations after large disturbances are crucial for 
understanding the susceptibility of ecosystems to 
extreme events such as those expected under altered 
climate scenarios (for example, Katz and Brown 
(1992) and works cited). 



Lake microbial communities are resilient 

A Shade et al 



2154 

Microbial communities are at the heart of all 
ecosystem functions, and thus their responses to 
disturbances may influence ecosystem recovery 
(Allison and Martiny, 2008). Some have advised that 
microbial variables be included in predictive models 
of ecosystem change (for example, Arhonditsis and 
Brett, 2004; McGuire and Treseder, 2010; Sarmento 
et al., 2010). But, before this can occur, it is necessary 
to understand the robustness of microbial assem- 
blages to disturbances (Allison and Martiny, 2008). 
In ecology, community robustness is comprised of 
resistance, defined as the ability to withstand change 
in the face of a disturbance, and resilience, the 
pace of recovery (if any) following a disturbance 
(Pimm, 1984). 

Water column mixing is a disturbance to microbial 
communities because it disrupts the physical- 
chemical gradients created by thermal stratification 
known to define niches for microorganisms (for 
example, Heaney and Tailing, 1980; Vincent et al., 
1984; Ovreas et al., 1997; Cytryn et al., 2000; 
Fenchel and Finlay, 2008). From an ecological 
perspective, bacterioplankton responses to mixing 
may provide unique insight into understanding 
disturbances. As it disrupts gradients in the water 
column and occurs seasonally, mixing affects spatial 
and temporal drivers of microbial communities. 
Further, both microbial and environmental res- 
ponses to mixing are tractable. High-throughput 
sequencing and fingerprinting methods can be used 
to assess changes in microbial community diversity, 
while environmental sensors can be used to quantify 
spatial and temporal changes in the environment. 
The alignment of these measurements allows a 
comprehensive perspective of microbial responses 
to an ecosystem-level disturbance. Therefore, lake 
mixing and microbial communities together create 
an interesting model system for understanding 
disturbance ecology. 

The work described here builds on a series of 
studies focused on the response of freshwater 
microbial communities to mixing (Shade et al., 
2007, 2008, 2010a, b, 2011; Jones et al., 2008). 
In an earlier study, we observed that patterns of 
bacterial community succession in a temperate 
eutrophic lake were linked to spring and fall mixing 
events (Shade et al., 2007), suggesting an impor- 
tance of lake mixing for community dynamics. In a 
second study of a small, darkly stained sub-tropical 
lake, we observed a surprising recovery among 
bacterial communities following mixing events 
caused by typhoons, with repeatable trajectories 
of community change across typhoon events and 
across years (Jones et al., 2008; Shade et al., 2010b). 
The results of these two surveys prompted us to 
ask whether bacteria in temperate lakes that mix 
seasonally would be similarly resilient to an 
episodic mixing event. Therefore, we conducted an 
enclosure experiment in a temperate humic lake to 
evaluate the effect of mixing and specific mixing- 
associated environmental changes on communities in 



the epilimnion (near surface stratum) and hypolim- 
nion (bottom stratum) (Shade et al., 2011). Both 
epilimnion and hypolimnion communities returned 
to their control composition within 10 days, demon- 
strating community robustness despite differences in 
initial compositions and in disturbance characteris- 
tics. Finally, we conducted a reciprocal transplant 
experiment to ask how epilimnion and hypolimnion 
bacterial communities responded to environmental 
conditions in the opposite strata (Shade et al., 2010a). 
Consequently, we identified possible 'generalist' 
community members, persistent in both epilimnion 
and hypolimnion, that could serve as pioneer species 
after mixing. The results from the transplant experi- 
ment suggested that lake bacterial communities 
harbor members that initiate post-disturbance succes- 
sion and eventual recovery. 

Following from the results of this series of studies, 
we asked whether freshwater bacterial communities 
would be robust (resistant and resilient) to an 
artificial mixing event that occurred during summer 
stratification, which is a time when communities are 
compositionally distinct from communities that 
occur during seasonal spring and fall mixing events. 
Mixing at the peak of summer stratification is a 
disturbance that would not happen naturally, and 
therefore represents an extreme event for summer 
microbial communities. Thus, this disturbance 
experiment provided a unique opportunity to 
measure both community and ecosystem responses 
to a large, unlikely disturbance in a relatively well- 
studied system. Summer mixing also allowed us to 
observe changes in the bacterial community separate 
from other seasonal changes that are concurrent with 
spring and autumn mixing (for example, ice-off). 
Further, a whole-ecosystem experiment avoided 
limitations of mesocosm and enclosure experiments 
(Schindler, 1998). One common criticism of meso- 
cosms is that they fail to reproduce conditions in 
nature and provide limited insight, and so it was 
possible that the community responses we previously 
observed in mesocosm experiments might not occur 
in the whole-lake. For these reasons, we wanted to 
pursue a whole-ecosystem experiment to understand 
microbial community responses to an unlikely, 
extreme mixing event. 

Here we report the bacterial community responses 
to the disturbance of an entire small temperate lake, 
North Sparkling Bog (Boulder Junction, WI, USA, 
Supplementary Table 1). We assessed bacterioplankton 
communities using 16S ribosomal RNA gene tag 
pyrosequencing to gain an unprecedented depth of 
sampling of the lake microbial diversity. 



Materials and Methods 

Experiment, site, environmental sensors and sampling 
The study lake, North Sparkling Bog (46°00'17.47' N, 
89°42'18.37' W), is darkly stained with high dissolved 
organic carbon (average 9.5 mgl" 1 ). Though this is 
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a relatively shallow lake (4.5 m maximum depth), 
the darkly stained color and surface area (0.46 ha) to 
depth ratio (1022 m) causes strong stratifica- 
tion during most of the ice-free season, which 
creates strong vertical gradients in temperature 
and dissolved oxygen, and overall high stability. 
North Sparkling Bog is a Global Lakes Ecological 
Observatory Network site. It has a buoy instru- 
mented with environmental sensors (see below); 
sensor data are available online (www.gleon.org) 
and lake characteristics are given in Supplemental 
Table 1. Whole-lake destratification was performed 
by mixing epilimnion and hypolimnion waters 
with the aid of buoyant forces (Read et al., 2011). 
Briefly, two rubber membranes with inflatable 
tubing on the outer rim (8.25 m diameter) were 
oscillated in the water column. When the tubing was 
filled with compressed air, the membranes rose to 
the lake surface. After air was released from the 
tubing, the membrane sank to the bottom. This 
repeated process of inflation and deflation caused 
turbulence and internal waves, and eventually 
resulted in complete mixing. 

The lake was stratified in oxygen and temperature 
before the artificial mixing that began on 02 July 

2008. An autoprofiling sonde (Yellow Springs 
Instruments, Yellow Springs, OH, USA; 6600 V2-4) 
measured hourly profiles of dissolved oxygen and 
temperature, and these observations were used to 
assess progression of lake mixing (data not shown). 
The buoy was also outfitted with a thermistor string 
(Apprise Technologies, Inc., Duluth, MN, USA) and 
a sensor that measured aqueous and atmospheric 
carbon dioxide concentrations every 2h (Browne, 
2004). The mixing manipulation continued until the 
surface and bottom-water temperatures differed by 
less than the accuracy of the thermistors (0.2 °C), on 
10 July 2008 (the last day of mixing, referred to 
throughout as Mix), signifying that complete mixing 
had been achieved. 

Samples were collected weekly or twice weekly 
during June, July and August of 2007, 2008 and 

2009, respectively, for standard water chemistry and 
microbial community analyses. Manual measure- 
ments of oxygen and temperature were collected at 
every sampling time (Yellow Springs Instrument 
550 A) and used to identify the thermocline. Inte- 
grated epilimnion and hypolimnion samples were 
collected as previously described (Shade et al., 
2008). In 2008, integrated thermal-layer samples 
were collected more frequently (usually daily) 
during and after the whole-lake manipulation. In 
2008 during the experimental period, we also 
collected depth-discrete epilimnion (0 m) and hypo- 
limnion (4 m) water samples using a peristaltic 
pump that was rinsed between each sample. 
Depth-discrete samples were collected in this 
manner even when the lake was completely mixed 
(and therefore no epilimnion or hypolimnion layers 
existed). Thus, on the day that complete mixing was 
achieved (Mix), we considered the 0m and 4 m 
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water samples to represent epilimnion and hypo- 
limnion water, respectively, even though the tem- 
perature throughout the water column was 
homogeneous. Collected water was preserved for 
water chemistry (total and dissolved nitrogen, total 
and dissolved phosphorus, dissolved organic and 
inorganic carbon, water color) using the North 
Temperate Lakes Long Term Ecological Research 
standard protocols (http://lter.limnology.wisc.edu/). 
Water for bacterial community analysis was filtered 
onto a 0.2 |im polyethersulfone membrane (Pall, 
New York, NY, USA), frozen at -20°C at the field 
station, and then transferred to — 80 °C at the main 
laboratory in Madison. 

During 2008, we also measured concentrations of 
several terminal electron acceptors and donors that 
are hallmarks of distinct biogeochemical transfor- 
mations mediated by microorganisms under anoxic 
conditions (for example, iron reduction, sulfate 
reduction and methanogenesis), throughout the 
ice-free season. Data presented here were collected 
on 01 July (day - 9), 10 July (Mix), 13 July (day 3), 
15 July (day 5), 17 July (day 7), 19 July (day 9), 21 
July (day 11) and 30 July (day 20). All parameters 
were measured at 4 m. We also measured five-point 
water column profiles of carbon dioxide and 
methane. For the carbon dioxide and methane 
profiles, one sample was always collected at 
0 m, one sample was always collected at 4 m to 
match the microbial samples, and three additional 
samples were collected at intermediate depths 
that best characterized the inflection points in 
dissolved oxygen over the oxycline. These middle 
three sample locations varied with the thermal 
and oxygen structure of the water column, and 
therefore were not the same depths across all 
sampling times. When the water column was mixed 
(no oxygen inflection point was observed), 
methane and carbon dioxide samples were collected 
equidistant from the surface and bottom at 0, 1, 2, 3 
and 4 m. 



Water chemistry 

Water chemistry analyses were performed by the 
North Temperate Lakes Long Term Ecological 
Research, according to their standard protocols 
(http://lter.limnology.wisc.edu/), except for iron, 
sulfide, carbon dioxide and methane. Samples for 
dissolved Fe were acidified with HC1 to pH<2 and 
analyzed colorimetrically using ferrozine (Stookey, 
1970). Samples for dissolved sulfide were stabilized 
by addition of 1% (wtvol -1 ) zinc acetate (to convert 
dissolved sulfide to ZnS), and analyzed colorime- 
trically via the methylene blue method (Cline, 1969). 
Carbon dioxide and methane were measured using 
gas chromatography (Shimadzu GC-14A, Kyoto, 
Japan), equipped with a flame ionization detector 
and a methanizer, which allowed for analysis of 
both CH 4 and C0 2 (which is converted to CH 4 by a 
methanizer) in a single injection. 
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Bacterial community sequencing and fingerprinting 
Total DNA was extracted from frozen filters as 
described using the FastPrep kit (Qbiogene, Carlsbad, 
CA, USA) with minor modifications (Fisher and 
Triplett, 1999). To track bacterial community 
responses to mixing, we employed a bar-coded 
pyrosequencing approach targeting the V2 variable 
region of the 16S ribosomal RNA gene (Fierer et al., 
2008; Hamady et al., 2008) to process samples 
collected from the epilimnion (0 m) and hypolimnion 
(4 m) on 01 July (day - 9), 10 July (Mix), 13 July (day 
3), 15 July (day 5), 17 July (day 7), 19 July (day 9), 21 
July (day 11) and 30 July (day 20). Pyrosequencing 
was performed using previously-described protocols 
for sequencing (Fierer et al., 2008), quality control, 
phylogenetic assignment (Hamady et al., 2010) and 
calculation of phylogenetic diversity (Faith, 1992). 
The sequencing targeted 250 base-pairs of the 16S 
ribosomal RNA gene using the primers Pyro 27f and 
338R. We recovered an average of 1441 sequences per 
sample (range 1248-1672). The phylogenetic tree was 
built using FastTree (Price et al., 2009). Analyses 
were performed with QIIME (Caporaso et al., 2010), 
using default parameters including the RDP classifier 
(Wang et al., 2007) for taxonomy assignment and cd- 
hit (Li and Godzik, 2006) with a 97% threshold for 
operational taxonomic unit (OTU) picking. Commu- 
nities were standardized by total number of reads. 

Rarefied Faith's phlyogenetic diversity (PD, 
(Faith, 1992)) and rarefied richness were calculated 
in QIIME. Pielou's Evenness, /, was calculated based 
on Shannon diversity, H [J=H/\og (richness)). The 
Grubb's outlier test (outliers package (Komsta, 
2011)) in the R environment for statistical comput- 
ing (R Development Core Team, 2010) was used to 
determine if the day of complete mixing values of 
richness, Pielou's evenness or PD were significantly 
different from the rest of the observations in the 
series. Grubb's test calculates a score for the outlier 
by subtracting the proposed outlier value from the 
set mean and then dividing by the s.d. [G statistic), 
and also by comparing the ratio of variances of the 
set with and without the proposed outlier included 
[U statistic). A paired, two-sided Mest was used to 
determine if epilimnion and hypolimnion trajec- 
tories of richness, Pielou's evenness or PD were 
significantly different from each other. 

Non-metric multidimensional scaling analysis 
was performed with UniFrac weighted distances 
(Lozupone et al., 2007) vegan package in R (Oksanen 
et al., 2011). The results are not sensitive to choice 
of metric: in this case, weighted UniFrac performed 
similarly to unweighted UniFrac and also Bray- 
Curtis similarity, suggesting that for this ecosystem 
patterns in community membership and community 
structure are highly correlated (see Supplemental 
Online Materials). The R environment was used for 
calculation of evenness and creation of bar charts 
and heatmaps by phyla. The vegan package in R was 
used for correspondence analysis, Mantel tests and 
Procrustes analysis. 



For the evaluation of the species abundance 
distributions (Magurran, 2007), we fitted log-normal 
models to the ranked abundance of each taxon in the 
pyrosequencing data, assuming an underlying Pois- 
son distribution (vegan package in R). Null, pre- 
emptive, log-normal, Zipf and Zipf-Mandelbrot 
models were each fit to the data (Wilson, 1991). 
Null and pre-emptive models did not produce good 
fits, as evaluated by Akaike information criterion, 
Bayesian information criterion and deviance. How- 
ever, the remaining models had comparable fits. The 
log-normal model was selected because it made the 
fewest assumptions about the underlying ecology of 
the data. Species abundance distributions for other 
models reproduced the same shift in the community 
structure at the day of mixing. 

To examine bacterial community change before, 
during and after the mixing event, we used 
automated ribosomal intergenic spacer analysis 
(ARISA) to fingerprint bacterial communities from 
integrated epilimnion and integrated hypolimnion 
water samples (Fisher and Triplett (1999) collected 
during June, July and August, with modifications 
provided in Shade et al. (2011)). This resulted in 
ARISA profiles for the following samples: 17 in 
2007, 60 samples in 2008 and 17 in 2009. Note that 
this integrated sampling strategy was different than 
what was used for the depth-discrete samples that 
were pyrosequenced. The demarcation of epilim- 
nion and hypolimnion samples was determined by 
thermal profiles measured at the time of sampling, 
and the thermocline was included in the epilimnion 
sample. Thus, integrated epilimnion and integrated 
hypolimnion together included the whole-lake 
column. 

ARISA profiles from integrated thermal-layer 
samples were used to determine to what extent 
physical displacement of OTUs via water motion 
could explain observed changes in community 
composition, during the July 2008 manipulation. 
We compared day — 10 samples (30 June 2008) with 
Mix (10 July 2008) to assess the proportion of Mix 
community members that could be attributed to 
OTUs present in the epilimnion or hypolimnion 
before the water column was mixed. Any commu- 
nity members that were present in the Mix sample 
but were not observed in the pre-mixing samples 
were designated as 'Unknown origin.' We used a 
Venn analysis, pie charts and hierarchical clustering 
with a heatmap of OTU occurrences to visualize 
differences in pre-mix and Mix communities (R 
Development Core Team, 2010). 

Next, we performed simulations of community 
richness and evenness to determine the probability 
of observing the Mix community given the pre-mix 
composition. For the simulation, we used re-sam- 
pling without replacement of integrated epilimnion 
and integrated hypolimnion communities in 'mix- 
ing scenarios' of the following proportions (epilim- 
nion/hypolimnion contribution): 10/90, 20/80, 
30/70, 40/60, 50/50, 60/40, 70/30, 80/20 and 90/10 
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(Supplemental Methods). For 1000 resamplings at 
each scenario, we report a P-value to describe how 
often we encountered the observed Mix richness 
and evenness in the simulated re-sampling. Using 
the ARISA data from the integrated water samples, 
we could evaluate whether the community compo- 
sition at Mix could be explained by different degrees 
of physical displacement of community members. 
Richness was compared rounded with the nearest 
whole number, and evenness was compared with 
three significant digits. More details for this analysis 
are provided in Supplemental Materials. 

We also used ARISA fingerprints to compare inter- 
annual dynamics in North Sparkling Bog the year 
before (2007, June through October) and the year 
after the manipulation (2009, June and July). These 
samples were fingerprinted from integrated epilim- 
nion and integrated hypolimnion water. If 
the lake bacterial communities recovered after the 
manipulation, we would expect the dynamics after 
the manipulation to be generally consistent across 
years because of the seasonality previously observed 
in North Temperate bog bacterial communities (Kent 
et al., 2007; Shade et al., 2008). To assess commu- 
nity variability, we calculated permuted analysis of 
betadispersion of pairwise Bray-Curtis similarities 
(vegan package (Oksanen et al., 2011) and 
(Anderson et al., 2006)). This analysis calculates 
dispersion from a mean compositional centroid 
within a categorical group. In this case, groups were 
months within years (for example, July 2007). As 
epilimnion and hypolimnion communities vary 
differently through time (Shade et al., 2007), 
epilimnion and hypolimnion samples were ana- 
lyzed separately. 

ARISA fingerprints were further used to calculate 
a rate of compositional change [a.k.a. similarity- 
decay) in the community. The slope of the log-linear 
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regression between Bray-Curtis dissimilarity and 
days between observations is a rate of compositional 
change (Nekola and White, 1999; Soininen et al., 
2007); this is similar to the strategy used to assess 
rate of change in the previous mesocosm experiment 
(Shade et al., 2011). We compared rates across the 
2007-2009 June-August time series to the 2008 
post-mixing time series (July-August), and also to 
the rate of change in temperature over the same time 
period. The temperature rate was assessed in the 
same way, by regressing the change in temperature 
against the time between observations. 

Results 

Physical and chemical changes after mixing 
Artificial mixing continued for 8 days until tem- 
perature was the same throughout the water column. 
As compared with the previous ice-free season 
(2007), mixing disrupted the normal physical and 
chemical structure of the lake, creating a pulse 
disturbance and exposing bacterial communities to 
markedly different environmental conditions. These 
different conditions included increased bottom- 
water temperature from 6 to 20°C and increased 
dissolved oxygen concentration from below detec- 
tion ( < 0.01 mg 1 " a ) to near 3 mg 1 ~ 1 (Figures la and b). 
The vertical oxygen gradient was homogenized 
gradually, but then returned within 2 days to its 
pre-manipulation anoxic conditions (12 July 2008, 
day 2) (Figures lc, d). 

We consider 'recovery' as the point in which the 
biology or chemistry matched the pre-mixing time 
point, day - 9. During the recovery period of the 
lake, stratification was relatively weak and more 
susceptible to weather-induced mixing, which is not 
typical for this typically strongly stratified lake 
(Figure la). We observed a spontaneous secondary 




Jun Jul Aug Sep Oct Nov Jun Jul Aug Sep Oct Nov 



Figure 1 Temperature and dissolved oxygen profiles of North Sparkling Bog during the year before mixing (2007; a, c) and the year of 
the mixing manipulation (2008; b, d). Profiles were observed with a hand-held dissolved oxygen and temperature sensor. Mixing began 
on 02 July 2008 (day — 8) and ended on 10 July 2008 (day 0, Mix) when homogenous temperature was observed. Black ticks above panels 
indicate sampling points. 
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mixing event 4 days after treatment had ceased 
(day 4). This secondary event occurred because 
overall the manipulation warmed the lake and 
created nearly isothermal conditions. The lack of 
thermal stability in the water column made the lake 
unusually susceptible to weather events that encou- 
rage mixing, such as temperature changes or wind. 
To the best of our knowledge (and based on 
inspection of thermal profiles collected from the 
time that instrumented buoys were first deployed on 
similar lakes, available at http://www.gleon.org), a 
spontaneous, full water column mixing event has 
never been observed in mid-summer in North 
Sparkling Bog. 

To understand any possible functional microbial 
responses that may be attributed to lake mixing, we 
measured the concentrations of various hypolim- 
nion respiration products over the mixing experi- 
ment. Measurements of hypolimnion iron (Fe 2+ and 
Fe Total ), sulfur (H 2 S) and nitrogen (NHJ respiration 
products decreased to a minimum on the day 
complete mixing was achieved (Mix), and then 
rapidly increased, generally following the overall 
response of dissolved oxygen (Figure 2a). S0 4 had a 
Tag' response that correlated with the post-mixing 
increase in H 2 S. Additionally, Fe 3+ and N0 3 did not 
vary substantially over the course of the experiment, 
and were both present in comparably low 
concentrations. 

We also considered the dynamics of gases C0 2 and 
CH 4 , two microbial products, over the course of the 



mixing experiment. C0 2 concentrations in the sur- 
face waters during the mixing experiment were 
markedly higher than C0 2 concentrations in the 
month prior and were more comparable with 
concentrations observed in the autumn during fall 
mixing (Figure 2b). Methane and C0 2 , two products 
of microbial respiration that often accumulate in the 
hypolimnia of stratified lakes, were homogenous 
throughout the water column at low concentrations 
at Mix, but re-accumulated at the lowest depth (4 m) 
by day 20 post-mixing (Figures 2c-e). 



Comparing patterns in ARISA fingerprinting to 
patterns in pyrosequencing 

We first compared bacterial community composition 
measurements made with automated ribosomal 
intergenic spacer analysis (ARISA) to those gener- 
ated via tag pyrosequencing from the same depth- 
discrete Om and 4 m July 2008 water samples, to 
ensure that they revealed the same community 
patterns. Using Bray-Curtis similarity for the ARISA 
fingerprints and weighted UniFrac distance for the 
pyrosequencing data, we observed similar patterns 
and results across both data sets (Mantel Pearson's 
i? = 0.79, P= 0.001 on 1000 permutations). We found 
a similar agreement across data sets by Procrustes 
correlation between correspondence analyses sam- 
ple scores [R = 0.75, P= 0.003 on 1000 permuta- 
tions). This analysis implied that community 
patterns generally were consistent across these two 
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Figure 2 Changes in measured water chemistry before and after lake mixing, (a) Microbial oxidation-reduction chemistry measured at 
4 m depth (hypolimnion) over the mixing experiment. Pre-mixing, Mix and post-mixing time points are on the x axis, (b) Dissolved gas 
concentrations over the mixing experiment. Daily average carbon dioxide concentration in the surface waters of North Sparkling Bog 
during the ice-free period of 2008, measured by a carbon dioxide sensor deployed on the lake, (c) Methane (open symbols) and carbon 
dioxide (closed symbols) profiles for day — 9 (pre-manipulation), for (d) Mix (when complete mixing was achieved), and for (e) day 20 
(post-mixing), c-e were measured using standard analyses of depth-discrete, manually-collected samples. 
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methods, and that results from one method are 
relevant to understanding results from the other 
(Supplemental Materials). 



Changes in bacterial diversity after mixing 
To characterize community robustness, we analyzed 
important ecological aspects of community diver- 
sity: richness (number of distinct taxa), evenness 
(equitability of taxa), structure (representation of 
dominant and rare taxa), PD, the number and 
relationship of phylogenetic lineages) and composi- 
tion (relative abundance of specific taxa defined as 
phylotypes at the 97% sequence similarity level), 
using the tag pyrosequence data. Changes in com- 
munity composition, loss of species richness and 
gain in species richness relative to the pre-mix 
conditions (possible bloom or growth of popula- 
tions) were used as indicators of community 
responses to the disturbance. July 2008 (the month 
of the mixing experiment) samples from water 
collected at depth 0 m (epilimnion) and depth 4 m 
(hypolimnion) were used to assess immediate com- 
munity responses to mixing. 

We predicted that the epilimnion (0 m) and 
hypolimnion (4 m) communities would be indistin- 
guishable under mixed conditions. However, rich- 
ness and evenness were different, whereas PD was 
comparable across depths (Figure 3). This suggested 
that despite thermal homogeneity of the lake, each 
stratum maintained a separate community structure, 
possibly driven by species sorting by environmental 
drivers other than temperature. Such a situation 
could arise if the time scale of vertical water 
movement (likely driven by convective mixing 
rather than wind, (Read et ah, 2012)) was slower 
than growth or loss rates. 

For each aspect of diversity evaluated and for each 
depth (0 and 4 m) of the lake, the bacterial commu- 
nities recovered within days to pre-perturbation 
states (Figure 3, Figure 4). There was an opposite 
response in diversity in each of the two stratum 
(Figure 3). The epilimnion communities became 
richer and included taxa that were more phylogen- 
etically divergent after the mixing (Mix). All 
measures of epilimnion diversity peaked on day 3 
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Figure 4 Log-normal models fit to ranked taxa abundance to 
assess changes in the rarity and dominance of members in the (a) 
epilimnion (0m) and (b) hypolimnion (4 m). Red arrows below 
the x axis show the shift in community structure after mixing, and 
black arrows show recovery. Each line color is a community 
structure at a different sampling time. In the epilimnion, 
communities experienced a second, greater shift in structure on 
day 3, shown by the two red arrows. 
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Figure 3 Bacterial response to lake mixing, assessed with bar-coded pyrosequencing of the 16S rrn gene V1-V2 region and analyzed for 
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post-mixing, and then returned to pre-mixing levels 
(Figure 3, open symbols). Pielou's evenness on day 3 
post-mixing was detected as significantly different 
from the remaining July 2008 time series by the 
Grubb's outlier test (G = 2.04, [7=0.002, P< 0.001). 
Together, these data suggest that the maximal 
change of the epilimnion community to the mixing 
experiment occurred not on the day of complete 
mixing, but after a 3 -day lag. 

In comparison, the hypolimnion communities 
were less rich (Grubb's G=1.97, [7=0.072, 
P=0.006), contained fewer phylotypes (G=1.93, 
[7=0.102, P= 0.012), and were less even [G= 3.145, 
[7=0.017, P< 0.001) on Mix, but rapidly recovered 
to pre-mixing levels by day 3 (Figure 3, closed 
symbols). These analyses suggest that the hypolim- 
nion community changed maximally around Mix. 
As an overall comparison of the diversity time 
series, epilimnion and hypolimnion communities 
were different in evenness and PD, but they were not 
different in richness, suggesting that though the 
number of OTUs observed at either depth were 
comparable, each depth contained a different repre- 
sentation of taxa (evenness T=2.93 P<0.05; PD 
T= -3.64, P<0.05; richness P = 0.27) 

Species abundance distributions were used to 
investigate changes in the bacterial communities' 
structures (Magurran, 2007; McGill et al, 2007). 



There was a strong shift in the structures of 
both epilimnion and hypolimnion communities 
(Figure 4). The epilimnion became more even (taxa 
shared similar relative abundances) during and 
immediately after the manipulation and the largest 
shift was between Mix and day 3 (Figure 4a). The 
hypolimnion became less even (loss of richness and 
dominance by a few taxa, Figure 4b) during the 
manipulation, and largely recovered by day 3. The 
recovery of the shape of the species abundance 
distributions further supports community recovery. 

Changes in bacterial community composition after 
mixing 

We next assessed changes in epilimnion and 
hypolimnion composition (defined as the relative 
abundances of OTUs) before and after lake mixing, 
using tag pyrosequence data (Figure 5). The epilim- 
nion community experienced a large change 
between day — 9 pre-mixing and Mix, a second 
smaller change between Mix and day 3 post mixing, 
and then another large change between days 3 and 7 
post-mixing (Figure 5a). After day 7, the post-mixing 
epilimnion community was most similar in compo- 
sition to the pre-mixing state, suggesting epilimnion 
community recovery (weighted UniFrac distance 
between epilimnion day —9 and day 7 = 0.13, 
Supplementary Figure la). 




Day of experiment 



Figure 5 Changes in bacterial community composition assessed with pyrosequencing and visualized by non-metric multidimensional 
scaling analysis (NMDS) of weighted UniFrac similarities in the (a) epilimnion (0 m), (b) hypolimnion (4m), (c) epilimnion (red lines) 
and hypolimnion (blue lines). The communities change along trajectories shown by the dashed lines. Orange circles identify the pre-mix 
time point (day — 9), green circles identify the day complete mixing was achieved (Mix), and yellow circles identify the recovery time 
point (day 20). (d) Differences in weighted UniFrac distance between the epilimnion (0m) and hypolimnion (4m) through the 
experiment. Note that epilimnion and hypolimnion are operationally defined on Mix because strata do not exist when water column 
temperature is homogeneous. 
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The hypolimnion community in North Sparkling 
Bog changed directionally until it returned to a 
composition similar to pre-mixing by day 11 post- 
mixing (Figure 5b; weighted UniFrac distance 
between hypolimnion day —9 and day 11 = 0.17; 
Supplementary Figure lb). As described above, the 
epilimnion and hypolimnion communities became 
less distinct during mixing (weighted UniFrac 
distance = 0.15), but then diverged quickly as thermal 
stratification was re-established (Figures 5c, d). 
Epilimnion and hypolimnion communities remained 
15% different by weighted UniFrac distance, though 
the dissimilarity between communities at the 
two depths reached a minimum at Mix (Figure 5d). 
Further, similarity between the epilimnion day 



20 and day — 9 communities was greater than the 
similarity between the same days in the hypolimnion, 
suggesting that the epilimnion community recovered 
more completely than the hypolimnion. 

There were distinct responses to mixing among 
phylogenetically coherent groups defined at the 
phylum or Proteobacteria class level (Figure 6). 
In the epilimnion, Gamma-Proteobacteria bloomed 
briefly on day 3 post-mixing, possibly responding to 
nitrogen or other nutrients up-welled from the 
hypolimnion after a lag for growth (Supplementary 
Figure 2), suggesting that there were community 
changes that could be attributed to biology. In the 
hypolimnion, populations of Acidobacteria and 
Delta-Proteobacteria were below detection at Mix, 
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Figure 6 Phylum-level (and Proteobacterial class level) changes in epilimnion (0 m) and hypolimnion (4 m) composition through the 
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Project. Asterisks indicate that no members associated with the phylum were detected. In the x axis, orange circles identify the pre-mix 
time point (day —9), green circles identify the day complete mixing was achieved (Mix) and yellow circles identify the recovery time 
point (day 20). See Supplementary Figure 3 for a different perspective of phylum-level dynamics, normalized by total occurrences within 
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but re-appeared by day 3. Several members of these 
latter phyla are strict anaerobes, and may have either 
died or were protected within a micro-niche when 
the hypolimnion was oxygenated. Bacteriodetes also 
decreased in the hypolimnion and gradually recov- 
ered within 20 days. Other phyla, such as hypolimnion 
Chlorobi and epilimnion Firmicutes, had dynamics 
that did not obviously correspond to the mixing 
experiment. 

To compare similar patterns of change across 
phyla, we standardized each phylum's dynamics to 
its maximum observed abundance (for example, a 
dominant phylum that changed twofold would be 
scaled the same as a rare phylum that also changed 
twofold, Supplementary Figure 3). In the epilim- 
nion, members of less-abundant phyla generally 
responded to mixing more dramatically than domi- 
nant phyla. In the hypolimnion, some rare phyla (for 
example, Gamma-Proteobacteria and Firmicutes) 
bloomed after the mix, while some abundant phyla 
(for example, Bacteroidetes and Acidobacteria) 
dropped below detection on Mix and recovered by 
day 3 post-mixing. 



Understanding community changes at mixing: 
displacement and ecology 

To determine the separate contributions of physical 
displacement of cells and biological processes (that 
is, active migration or growth of OTUs), we used 
ARISA fingerprints of integrated epilimnion and 
integrated hypolimnion water samples (Figure 7). 
ARISA OTUs that were not detected before the 
manipulation (day - 10, note that this is 1 day 
earlier that the day — 9 time point used for 
pyrosequencing) were unique to the 'Epi Mix' and 
'Hypo Mix' communities (Figure 7a). These unique 
Mix OTUs were of unknown origin, and represented 
11% of the epilimnion and 12% of the hypolimnion 
communities (Figure 7b). Further, the epilimnion 
mix and hypolimnion mix communities were more 
similar to each other than they were to either the 
epilimnion or the hypolimnion community before 
mixing (Figure 7c). 

Finally, we performed re-sampling of the day — 10 
communities in various scenarios of mixing to 
assess how likely it was to observe the Mix 
communities given the day — 10 composition. We 
re-sampled communities one thousand times, 
including epilimnion and hypolimnion proportions 
of 10/90, 20/80, 30/70, 40/60 and 50/50, and 
assessed richness (assessed by the number of unique 
ARISA fragment lengths) and Pielou's evenness for 
each re-sampling. For these tests (a total of nine 
mixing scenarios of different proportions of epilim- 
nion and hypolimnion contribution), we never 
observed the actual ARISA-based richness and 
evenness after 1000 resamplings at each mixing 
scenario (Supplementary Figure 4). Together, these 
results suggest that there were changes in composi- 
tion that could not be attributed only to physical 



displacement of cells, but that may be instead 
explained by ecological processes. 



Long- and short-term effects on bacterial community 
dynamics 

To evaluate whether the mixing experiment caused 
prolonged changes in bacterial community 
dynamics, we assessed seasonal variability using 
community ARISA fingerprints from 2007 through 
July 2009. Because bacterial communities in north 
temperate bog lakes were previously shown to 
change seasonally (Kent et al., 2007; Shade et al., 
2008) we used permuted analysis of betadispersion 
of Bray-Curtis similarities to address whether there 
were overall monthly differences in variability 
between the year of mixing and the other 2 years. 
We detected no differences in overall variability 
across months and years (epilimnion P=0.35 
on 999 permutations, hypolimnion P=0.10, 
Supplementary Figure 5). This suggests that the 
mixing experiment in July 2008 did not have 
prolonged or permanent effects on bacterial com- 
munity dynamics. 

We then asked whether the mixing manipulation 
influenced the bacterial community rate of change, 
by comparing rates observed in the epilimnia and 
hypolimnia from 2007, 2008 and 2009. We calcu- 
lated the rate of change (similarity-decay) by regres- 
sing the Bray-Curtis dissimilarity (calculated using 
ARISA profiles) against the time interval between 
samples, for all possible pairs of samples within 
each time series, and recording the slope of the 
regression (Supplementary Figure 6). The hypolim- 
nion communities sampled in 2008 post-mixing 
changed faster than any of the other communities. 
The bacterial community rate of change was corre- 
lated with the rate of change in temperature 
(r=0.94, P= 0.005, 77 = 6). Taken together, these 
findings suggest that the elevated temperature in 
the hypolimnion contributed to a faster rate of 
community change than is normally observed in 
the cold bottom waters of this bog lake. 



Discussion 

We employed a whole-ecosystem manipulation to 
ask whether microbial communities could recover 
after a large disturbance, and to ask how responses 
of microbial communities compared with changes 
in environmental conditions. Bacterial communities 
were not resistant to summer mixing, but were 
highly resilient. The communities in the epilimnion 
and hypolimnion both changed as a result of the 
disturbance, but recovered by days 7 and 11 post- 
disturbance, respectively. Though microbial com- 
munities in the epilimnion and hypolimnion had 
contrasting responses to the disturbance (for exam- 
ple, epilimnion increased in diversity and the 
hypolimnion decreased), they had similar patterns 
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of recovery. Although the manipulation seemed to 
increase the rate of community change in the hypo- 
limnion during summer 2008, the communities did 
not exhibit prolonged changes in dynamics after the 
mixing experiment, as assessed by comparing over- 
all variability in the year of mixing to ARISA 
fingerprints from 2007 through July 2009. 

The richness and evenness of the integrated 
epilimnion and integrated hypolimnion at Mix were 
not the same as any simulated mixing scenarios, 
suggesting that vertical transport of cells in the 
water column was not enough to explain differences 
between the two layers (between day — 10 and Mix). 
We also observed changes in individual phylotype 
and OTU relative abundances that could not be 
explained by cell transport in the water column. 
For example, at mixing, Alpha-Proteobacteria tags 
increased in relative abundance in the hypolimnion, 
while other tags previously abundant in the hypo- 
limnion fell below detection (for example, Acido- 
bacteria, Delta-Proteobacteria). Several members of 
these latter phyla are strict anaerobes, and were 
either protected in micro-niches or possibly died 
when the hypolimnion was oxygenated. Interest- 
ingly, the Beta-Proteobacteria decreased in the 
epilimnion as a proportion of total tags, but changed 
minimally in the hypolimnion. Thus, we infer that 
biological processes (for example, growth, death, 
revival from a dormant state) contributed to the 
observed temporal dynamics. For example, dormant 
microorganisms persisting in a sediment seed bank 
may have a role in community resilience after 
disturbance (Jones and Lennon, 2010), but further 
study is needed to determine the role of inactive 
members, as well as extremely rare members, to 
overall community recovery. 

Hypolimnion response and recovery 
Hypolimnion communities changed more drasti- 
cally than epilimnion in composition, richness and 
phylogenetic diversity. The recovery of hypolim- 
nion composition coincided with the depletion of 
dissolved oxygen and the re-accumulation of anae- 
robic respiration products such as methane, carbon 
dioxide and other reduced molecules (for example, 
ferrous iron, ammonium and sulfide). As would be 
expected, accumulation of sulfide coincided with a 
decrease in sulfate, several days after Mix. Taken 
together, these observations are consistent with the 
transitions from aerobic to anaerobic respiration 
that we would predict from first-principles of 
biogeochemistry (for example, electron-donor and 
-acceptor pairs). Thus, the microorganisms and their 
respiration products helped to create the environ- 
mental conditions that eventually harbored a recov- 
ered hypolimnion community. This demonstrates a 
fundamental feedback loop for the hypolimnion 
community: microorganisms help to create the 
strong oxidation-reduction gradients, and this 
gradient then dictates the observed community 



robustness to mixing by acting as an environmental 
filter that drives species sorting. 

In our previous mesocosm experiment, the hypo- 
limnion community also recovered as oxygen con- 
centration decreased, despite that temperature 
immediately recovered to pre-disturbance condi- 
tions (Shade et al., 2011). Notably, the hypolimnion 
temperature was approximately 6 °C during the 
mesocosm experiment (the same as the ambient 
hypolimnion temperature) and the community 
recovered after 10 days, while in the whole-lake 
manipulation the temperature was 20 °C (elevated as 
compared with the pre-disturbance 6°C) and the 
community recovered after 11 days. This indicates 
that oxygen concentration is more important than 
temperature in determining hypolimnion commu- 
nity composition. This does not exclude the 
possibility that other, unmeasured environmental 
changes also were important in driving the micro- 
bial community, but oxygen proved to be the most 
important among those that we measured. However, 
elevated temperatures appeared to accelerate the 
hypolimnion community changes after mixing. 
Recent estimates of bacteriovore grazing rates in a 
similar small bog lake suggested bacterial growth 
rates were on the order of 0.33 day even at 9°C 
(Hahn et al., 2012). This is consistent with the 
community composition in the recovering hypolim- 
nion turning over during a period of days. 



Summer lake mixing as a novel, large disturbance 
We have provided several lines of evidence that our 
mixing experiment represented both a novel and 
extreme disturbance to lake bacteria. The first is the 
drastic changes in temperature and oxygen, key 
drivers of epilimnion and hypolimnion aquatic 
bacterial communities, respectively (Shade et al., 
2011). The temperature in the hypolimnion 
increased from 6 to 20 °C. The historical ice-free 
hypolimnion temperature has likely never exceeded 
10 °C. The oxygen in the hypolimnion was increased 
from 0 to near 3 mg/1. Stable bog lakes such as North 
Sparkling Bog remain anaerobic in the bottom 
waters with the exception of seasonal mixing (fall 
and spring). Both of these environmental changes 
were major deviations from seasonal norms. Second, 
the known seasonality of the microbial communities 
in North Sparkling Bog (Kent et al., 2007) provides 
additional evidence for the novelty of the mixing 
event. Because North Sparkling Bog bacterial com- 
munities typically follow a temperate lake seasonal 
trajectory, we have observed that community com- 
position in July is very different from the composi- 
tion observed during spring and fall mixing. 
Therefore, the mixing disturbance was novel to 
summer 'July' communities. Finally, the duration of 
the mixing event presents an extreme disturbance 
for this ecosystem. After 8 days, complete mixing 
was achieved in the lake, whereas seasonal spring 
and autumn mixing occur after weeks to months of 
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gradual cooling of the epilimnion and periodic 
entrainment of these cooled waters into the hypo- 
limnion, until temperature homogeneity is reached 
through to the sediment of the lake. The experi- 
mental mixing occurred over a much shorter period 
than seasonal mixing. 

Conclusions 

Microorganisms as sentinels of change 
Our results have implications for understanding 
microbiological responses to environmental pertur- 
bations at both fine and coarse temporal scales. In 
this study, the microbial communities changed after 
the disturbance, exhibiting sensitivity to a pulse 
perturbation that occurred over a relatively short 
time scale (1 week). Further, the clear connection 
between biogeochemical conditions and community 
composition suggests that there are direct linkages 
between changes in environmental forcing, biogeo- 
chemical conditions, and microbial communities. 
The communities recovered within days, and the 
analysis of variability from 2007 to 2009 demon- 
strated that the monthly baseline in temporal 
community variation was relatively stable and may 
represent a deterministic point of return. Work by 
our group (Kent et al, 2007; Shade et al, 2007; 
Shade et ah, 2008) and others (Crump and 
Hobbie, 2005; Fuhrman et al, 2006; Nelson, 2008; 
Andersson et ah, 2010) has shown repeatable 
seasonal patterns in microbial communities inhabit- 
ing rivers, oceans, seas and lakes, providing cumu- 
lative evidence of a consistent microbial community 
'baseline' that spans many aquatic systems. In 
ecology, increases in ecosystem variance have been 
suggested as indicators of ecological regime shifts 
(for example, Carpenter and Brock, 2006). Therefore, 
gradual changes in aquatic microbial communities 
that exceed normal variability may correspond to 
press disturbances, such as those expected with 
gradual global change (for example, increases in 
mean temperature). Meanwhile, more dramatic 
microbial shifts may signal pulse disturbances, such 
as extremes in precipitation or temperature asso- 
ciated with expected increases in climate variability 
(Folland et ah, 2002). The current work and these 
previous studies together indicate that consistent 
time series of aquatic microbial communities may 
serve as robust sentinels over multiple time scales, 
from local pulse disturbances to global presses. 
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